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LONG-TERM  GOAL 

My  longterm  goal  is  to  develop  physics-based  models  to  predict  nearshore  morphology  changes, 
and  to  test  the  models  against  available  field  data. 

OBJECTIVES 

The  objectives  of  the  present  project  are  to: 

1.  Incorporate  a  sediment  transport  and  bed  morphology  capability  in  the  wave-induced  circulation 
model  of  Ozkan-Haller  and  Kirby  (1997,  1999) 

2.  Use  the  resulting  model  to  study  the  growth  to  finite  amplitude  of  bottom  perturbations  on  initially 
longshore-uniform  planar  or  barred  beaches 

3.  Use  the  model  to  investigate  the  evolution  of  three-dimensional  bed  features  under  specific  wave 
conditions,  in  comparison  to  ARGUS  video  results  for  the  Duck  FRF  site. 

4.  Begin  the  development  of  an  instantaneous  sediment  transport  capability  (on  a  wave  by  wave 
basis)  in  the  Boussinesq  model  of  Wei  et  al  (1995). 

APPROACH 

Our  approach  is  to  use  a  robust  numerical  code  for  the  modeling  of  2D  or  quasi-3D  wave-driven 
nearshore  circulation  as  the  basis  for  computing  the  local  wave-averaged  sediment  transport  rate  and  the 
resulting  evolution  of  nearshore  morphology.  The  model  results  will  be  limited  in  accuracy  by  both  the 
accuracy  of  input  wave  information  and  by  the  accuracy  of  the  local  sediment  transport  model.  By  using 
data  from  the  FRF  ARGUS  station  and  related  wave  information  from  in  situ  pressure  gage  arrays,  we 
seek  to  determine  whether  the  model  will  reproduce  qualitative  shifts  in  overall  bed  geometry  which 
have  been  observed  to  correlate  with  shifts  in  wave  conditions  (Fippman  and  Holman,  1990). 

WORK  COMPLETED 

We  have  incorporated  a  sediment  transport  formulation  and  bed  morphology  calculation  into  the 
nearshore  circulation  model  of  Ozkan-Haller  and  Kirby  (1997,  1999),  and  have  begun  to  use  the  model 
to  study  evolution  of  large  amplitude  bed  forms  in  river  environments  as  well  as  the  nearshore  environ¬ 
ments.  The  study  of  river  environments  has  proven  to  be  a  useful  intermediate  step,  as  it  has  allowed  us 
to  test  the  robustness  of  the  computational  model  before  concentrating  on  the  accuracy  of  the  sediment 
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transport  formulation  after  waves  are  included.  River  environments  provide  striking  examples  of  large 
amplitude  bed  features,  including  alternating  bar  configurations  with  strong  shock-like  structures.  We 
have  begun  an  evaluation  of  the  long-term  behavior  of  various  unstable  bedform  patterns,  whose  lin¬ 
ear  stability  analyses  have  been  described  in  the  literature  (Falques  et  al,  1996;  Schielen  et  al,  1993). 
Work  on  extending  the  model  to  incorporate  a  Bagnold-Bailard-Bowen  type  model  for  wave-induced 
transport  is  underway. 

RESULTS 

Results  to  date  have  concentrated  on  calculations  based  on  sediment  transport  modeling  for  slowly 
varying  unidirectional  flows,  as  pursued  in  Falques  et  al.  (1996)  for  coastal  environments,  or  by  Parker 
(1976),  Colombini  et  al  (1987),  Schielen  et  al  (1993)  and  others  for  river  environments.  In  particular, 
we  have  concentrated  recently  on  the  evolution  of  alternate  bars  in  river  environments.  The  growth  of 
these  bars  as  a  bed  instability  has  been  understood  as  a  linear  instability  process  since  Parker  (1976), 
and  Colombini  et  al  (1987)  showed  in  a  next-higher-order  correction  that  the  bar  crests  could  become 
oblique  to  the  channel  sides.  However,  it  has  long  been  understood  from  observations  that  alternate  bars 
can  grow  into  very  large  amplitude  structures  with  very  abrupt,  down-channel  facing  depth  increases, 
taking  on  a  shock-like  appearance.  Examples  in  Chang  et  al  (1971)  serve  as  a  guideline  for  testing.  An 
example  from  our  own  numerical  work  is  illustrated  in  Figures  1  and  2.  Figure  1  illustrates  the  evolution 
of  a  single  along-channel  wavelength  of  a  nearly-fastest  growing  linear  perturbation  in  a  channel  200 
m  wide  and  5  m  deep.  In  these  figures,  along-channel  is  oriented  up  the  page.  The  initially  sinusoidal 
bars  shown  at  t  =  0  take  on  an  oblique  character  and  evolve  into  shock-like  structures  fairly  rapidly, 
by  t  =  9 hr  as  illustrated  on  the  rightmost  panel.  Figure  2  illustrates  the  continuation  of  this  process 
at  much  later  time,  where  the  shock  structure  has  equilabrated  and  is  propagating  along  the  channel. 
Depths  in  the  scour  holes  immediately  downstream  of  the  bar  crests  can  be  twice  the  initial  unperturbed 
channel  depth. 

Features  such  as  these  show  up  in  a  wide  range  of  test  calculations  with  a  variety  of  sediment 
transport  formulations  and  channel  geometries,  and  thus  are  likely  to  self-organizing  features  to  some 
degree.  However,  the  quantitative  accuracy  of  the  morphology  predictions  here  is  almost  certainly 
dependent  on  the  chosen  sediment  transport  relation  as  well  as  the  limitation  to  depth-integrated  2-D 
hydrodynamics. 

The  model  has  subsequently  been  applied  to  coastal  beaches  in  two  scenarios.  In  the  first,  we  have 
used  a  formulation  similar  to  that  employed  by  Falques  et  al  (1996),  in  which  sediment  transport  is 
driven  primarily  by  the  mean  flow  and  wave-induced  effects  (such  as  cross-shore  transport  in  the  pres¬ 
ence  ofa  longshore  current)  are  absent.  Figures  3  and  4  show  the  early  and  later  stages  of  growth  of  an 
initial  perturbation  similar  in  form  to  the  predictions  of  the  linear  stability  theory  of  Falques  et  al.  Ini¬ 
tially,  the  growth  of  relatively  well  organized  oblique  bedforms  is  observed.  This  organized  evolution 
does  not  persist  in  the  present  simulations,  however,  where  the  bed  evolves  further  to  a  somewhat  cus¬ 
pate  from  with  a  hint  of  the  formation  of  rip  channels.  Further  interpretation  of  these  results  requires  the 
application  of  a  much  more  accurate  wave  model  which  accounts  for  wave-current  interaction  effects. 

Work  has  also  been  started  on  performing  similar  simulations  using  the  Bagnold-Bailard-Bowen 
transport  formulation  or  more  modern  variants  incorporating  acceleration  effects.  At  the  end  of  this 
project,  we  are  able  to  predict  the  formation  of  a  storm  bar  under  strong  erosional  conditions.  This 
work  is  being  continued  in  the  context  of  the  NOPP  Nearshore  Community  Model  project. 
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Figure  1:  Initial  growth  of  a  sinusoidal  alternating  bar  in  a  sloped  channel  The  initial  bar 
configuration  Is  based  on  linear  perturbation  analysis  and  has  an  along-channel  length  slightly 
longer  than  the  fastest  growing  mode.  Initial  channel  depth  is  5  meters ,  and  labelled  contours  are 

depths  in  meters. 
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Figure  2:  Later  stages  of  alternating  bar  evolution ,  showing  equilibration  and  along-slope  migration 
of  bars.  Note  deep  scour  holes  immediately  downstream  (or  up-page)  from  bar  crests,  as  well  as  the 
steadiness  of  the  zig-zagging  shock-like  structures  in  the  bed. 
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Figure  3:  Initial  growth  of  a  perturbed  bedform  corresponding  to  a  bed-surf  instability  described  by 
Falques  et  al  (1996).  Top  panels  shows  contours  of  bed  deformation  relative  to  a  planar  beach. 

Lower  panels  show  the  total  depth  profiles. 


Figure  4:  Later  stages  of  nearshore  bedform  evolution,  showing  formation  of  a  complex  but 

somewhat  cuspate  bedform. 
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IMPACT/APPLICATION 

This  work  is  intended  to  provide  a  process-based  model  for  predicting  nearshore  large  scale  bedform 
evolution  over  day  to  week  time  scales,  as  a  first  step  to  bridging  the  gap  between  short-time  physical 
predictions  and  longer  time  (annual  and  longer)  predictions  based  on  parameterized  physics.  This  gap 
is  still  immense,  and  will  only  be  closed  with  (1)  the  advent  of  both  faster  computers  and  with  (2)  the 
ability  to  assess  the  accuracy  of  local  sediment  transport  formulations  applied  over  reasonably  long 
time  scales.  We  are  mainly  concerned  with  the  second  aspect  of  this  problem. 

TRANSITIONS 

The  work  on  morphology  evolution  conducted  here  will  carry  over  directly  into  the  seabed  module 
work  in  the  NOPP  nearshore  project.  Tuba  Ozkan-Haller  is  working  very  closely  with  us,  with  model 
extensions  and  developments  being  shared  by  all  investigators. 

RELATED  PROJECTS 

(1)  N00014-99- 1-0490  ’’Prediction  of  the  Low  Frequency  Wave  Field  on  Open  Coastal  Beaches”, 
H.  Tuba  Ozkan-Haller,  Univ.  Mich.  Tuba  is  directly  involved  in  this  work,  and  will  be  in  principal 
charge  of  transitioning  results  on  the  morphology  model  here  to  NOPP  models.  We  are  also  using  her 
work  on  extensions  to  the  wave  driver  portion  of  the  hydrodynamic  code. 

(2)  N00014-99- 1-1051  (NOPP)  ’’Development  and  Verification  of  a  Comprehensive  Community 
Model  for  Physical  Processes  in  the  Nearshore  Ocean”,  James  T.  Kirby  et  al,  Univ.  of  Del.  The  work 
on  morphology  evolution  done  here  will  provide  a  foundation  for  work  in  this  area  in  the  NOPP  project. 
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